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A variational principle and conservation theorems 
in connexion with the generally relativistic 
Dirae equation* 


By B. E. Laurent 


Introduction 


In a recent paper Klein [1] has proposed a treatment of the gravitational field 
based on the use of the y and [’ matrices from the theory of the generally covariant 
Dirac equation (see Bargmann [2]). Klein introduces an action integral containing 


_ these matrices and the two sets are varied independently. This procedure bears 


a certain analogy to the Palatini method, in which the components of the metric 
tensor and of the affine connexion are varied independently. The Palatini method 
has proved to be well suited for attacking the problem of quantizing the gravita- 
tional field (Arnowitt and Deser [3]). One result obtained here is that the equations 
one obtains with Klein’s variation principle are the Einstein equations if the varia- 
tions of the y matrices are restricted so that 


hy pag 


where the g’ are multiples of the unit matrix, and are also the components of 
the tensor reciprocal to the metric tensor (see section [.4.). 

In section I the method with independent variations of the y and [’ matrices 
is extended to a gravitational field interacting with a Dirac field.1 (This corre- 
sponds to Weyl’s [4] extension of the Palatini method.) This is done in section L.5, 
where the new equations are obtained and compared with the Einstein and Dirac 
equations. The sections preceding this one are preparative in nature. Sections I.1-3 
constitute a summary of the Bargmann paper [2] with slight generalizations and 
section 1.4 is devoted to a comparative review of the usual variation principle 
with coupled variations of the y and the I’ matrices. 

In section II we derive an ‘‘energy-momentum complex” from the above-men- 
tioned Lagrangian without interaction. This is done in section II.1 following the 
procedure used by Moller in [5] (see also Laurent [6]); the resulting “complex” is 
very similar to Moller’s. In section II.2 the two are compared. Section II.3 deals 
briefly with a conservation identity connected with the S transformation. 


* The research upon which this paper is based was done under U.S. Air Force contract no. AF 
61(052)-47. “al 
1 ¥or simplicity the Dirac field for a zero-mass particle is chosen. 
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I. Field equations 


1.1. Notation 


In the generally relativistic theory of the Dirac equation one introduces the 
four 4x4 matrices y"(u=0, 1,2, 3), functions of the space-time coordinates, which 
together form a contravariant vector field. The anticommutators of these matrices 
are required to be multiples of the unit matrix,’ 


fy vy=2g""1, (1) 


and the tensor field g’” is identified with the reciprocal of the metric field 4,,. 
Indices are raised and lowered with the help of g’” and g,,. 
It is convenient to use the following definitions: 


se = sly): (2) 
_ aByd 7 3 
Vs= 542 VaYBY rs; (3) 
g= — Det (9), (4) 


where €*9 is the Levi-Civita symbol, antisymmetric in all indices and with e°”’ = 1. 


Under coordinate transformations s“” is a contravariant tensor, y; is a scalar density 
of weight 1, and g is a scalar density of weight 2. Equation (1) is covariant also 
under an arbitrary ‘‘S transformation’, 


(oh args, (5) 


where the matrix S generally is a function of the coordinates. 

All the quantities defined above transform in this manner under S transfor- 
mations. They are tensor-operators in Bargmann’s [2] terminology. In Appendix 
1 a few relations between these quantities are collected. 

The covariant derivative defined with help of the Ju Will, as usual, be denoted 
by a semicolon (e.g. y",,). Thus we have 


Jap:y = 0. (6) 
The y“,, and similar quantities derived from other tensor-operators are not gener- 


ally covariant under S transformations. We can, however, with help of a new 
matrix field, I’, construct a quantity 


UY ce a7 [y", P.] (7) 


(with corresponding equations for other tensor-operators) which is a new tensor- 
operator. For this purpose we let I, transform as a covariant vector under co- 
ordinate transformations and in the following manner under JS transformations: 


(T.)'= St? PeS- 877 Sys (8) 
* In the following, unit matrices will generally not be written out explicitly. 


264 


—o 


ARKIV FOR FYSIK. Bd 16 nr 25 


where the comma denotes the ordinary partial derivative. It is easy to check, with 
the help of (8), that y\ is, in fact, a tensor-operator and that the ordinary rule 
for differentiating a product applies also to the ‘“stroke-operation” on a product 


: of tensor-operators. We will give the operation (7) the name ‘extended covariant 
_ derivative’’ (e.c.d.). 


A new tensor-operator 
On 1h 6 Lae Lvl (9) 
can be constructed from a I field. It is easy to check directly, with the help 
of (8), that ®,, is, in fact, a tensor-operator. 
Finally, we notice that although I’, is not a tensor-operator, its variation, ole 
_ is. This can be used to express the variation of ®,, in the following convenient way! 
) ®,, 7 (6 Di oe (6 ins (10) 


which follows immediately from varying (9). 


1.2. The I” field 


In the following we shall use different I’, fields, all characterized by being co- 

variant vector fields and transforming according to (8) under S transformations. 

_ <Any [, field can be used for constructing an e.c.d. but, of course, one obtains 

a different e.c.d. for each I’, field. We shall distinguish between them by putting 
primes on the e.c.d. stroke; e.g. 


Ne th adhe iy} | eal (11) 
The field denoted by I) is defined through the equation 
yi» = 0. (12) 


This equation can be solved for the I, by expanding it in the 7 matrices of 
Appendix 1. The coefficients in this expansion may then be determined with the 
help of the relations (A1.9-12). The result is 


Ty = —2{4 Tr (Sap Vr) ¥* +3 Tr (Vp Var) 8 +g * Tr Ysa») Ys" — 
—497* Tr (ysyov'w) Ys} + %- (13) 


The quantity a, is here an arbitrary vector field transforming in the following 
manner under an S transformation 


(ay) = a, —}|S|-*| S|. (14) 
One can always, with the help of an S transformation, make? 


Vv ONY (15) 


where the C”,, are multiples of the unit matrix. It is often convenient to use 
this special system; then (13), for example, becomes 


1 This corresponds to Palatini’s formula giving the variation of the curvature tensor in terms 
of the variation of the affinity (see e.g. Moller [7] p. 334). 
2 ef. Bargmann’s [2] “distinguished”’ (ausgezeichnet) system. 
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has a (Ya Ypn)8"? + ay (16) 


if use is made of (A1. 8). Notice that Tr (y. ys.) is antisymmetric in « and 6 due 
to (1) and (6). ; 

A field ®;, is obtained from I’, through the general formula (9). It is easy 
to show, with the help of a locally inertial system and (15) (hence (16) may be 
used), that? 

Diy=} mie 8" + Ay — Buys (17) 
where R,,,, is the curvature tensor. It is immediately seen that this equation is 
covariant in all respects and so we can drop our restrictions. 

Using (A1.8) and (A1.10) we obtain from (17) 


Roouv = % TY (Sec ®,,). (18) 


The contracted form of this equation appears in Klein [1]. 


1.3. Dirac fields 


In Bargmann’s [2] formalism the Dirac fields are scalars transforming in the 
following manner under S transformations 


y=S*y. (19) 

We shall, with Bargmann, introduce a quantity « such that 
— yt = ay a? (20) 
and a*—a=0, (21) 


where the cross stands for hermitian conjugate. The matrix « is determined from 
(20) and (21) to within a real factor. The matrix « is taken to be a scalar, 
transforming under S transformations as follows 


ai = S* a 8. (22) 


One can easily check that (22) makes (20) and (21) covariant. A tensor with 
the property (22) is, in Bargmann’s terminology, a tensor-«, while a tensor with 
the property (19) is a tensor-y. 

We shall quite generally require from the e.c.d. that the ordinary rule for 
differentiation of a product applies and that it and hermitian conjugation are 
commuting operations. With the following definition for the e.c.d. of a tensor-« 


Ge Ont ely Ts x (23) 


the two rules check with (20) and (21) and a, is a new tensor-c. 
The relation 


op = 0 (24) 


does not contradict anything so far said about a. It gives a connexion between 
the undetermined fields in I’, and «. In the following we shall always assume 
(24) to be valid. 


* Write out the covariant derivative in (16) explicitly 


VBsv = VB.» — Ty Yo: 
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If the e.c.d. of a tensor-y is defined as follows 


, Pue=Vin — Lup = (25) 
- the two rules still apply, e.g. for products of the type yp’ «yy (which are tensor- 
operators). 


_- Using the two rules and (23), (24) and (25) we can further deduce the expression 
for the e.c.d. of the quantity p=y" a: 


Pu= Put Dy. (26) 

It is often convenient to use the quantities 
p=s~ (27) 
and P=" «a (28) 


as substitutions for y and @ (Klein [8]). 


1.4. The action containing T.. 


The action for the gravitational field coupled to a massless Dirac field is the 
following 


. [R42 Gy"qr,) dx, (29) 
where 
R=g' Rk 


and R is the curvature scalar: R= R**;,. 

To obtain the field equations from (29) one must vary g,,, p and g indepen- 
dently and set the corresponding variations of (29) equal to zero. Now, instead 
of varying the g,,, we can vary the y” in such a way that {y”,y’} is always a 
multiple of the unit matrix, i.e. 

Iya’, (30) 
where €”, is an infinitesimal mixed tensor field. The symmetric part of €", corre- 
sponds to a variation of g”,€”+€"=69", while the antisymmetric part corre- 
sponds to a class of § transformations and consequently gives no new field equations, 
since (29) is invariant under all S transformations. Hence we will in the following 
assume €”” to be symmetric. 

Using (18) and performing a partial integration we can replace (29) by 


[{ Tr 8 O29) — 2% (Gy? Qa Gray" p) } Ce. (31) 

If we introduce the explicit expressions for gy, and Gy, 
Pin = Pu Flan —Tue, (32) 
Pin = Gu 3 au G + GV, (33) 


into (31) we see that the middle terms in (32) and (33) cancel (Klein [8]). Con- 
sequently we need not consider them when performing the variations. 

The variation da, in (A2.4) is not determined by dy“, but we see that if we 
do not wish to obtain a vanishing current 
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Gy p=0 
as a part of the field equations, we must restrict the variation of I, so that da, =0, i.e. 
Tr (61) =0. (34) 


Equation (34) will be used throughout. 
The field equations are now obtained in a fairly straightforward manner. The 
variations of g and yield 


YG =0 (35) 
and Gruyi=0. - (36) 


In the first term of (31) the variation of ®,, gives no contribution due to (10) 
and (12), while the variation of 3%’ give us the following term in the variation of (31) 


— 4 f et, (M%,— hgh) d*e. 


Equations (A2.2) and (18) have been used here. 
From the second term of (31) we obtain the following term by varying y“ 


—2[(Gdy" Qu— Gnd y"p— Gly", él, } pl d*a. 


We notice here that the last term vanishes according to (A2.4) and (A 1.22). 
The use of (30) gives the final form to our expression, ie. 


— Hf [E (Gye Qvn— Fru VP + Gyn Prva -GraYug) dx 
and setting the variation of the total action (31) equal to zero gives us the Ein- 
stein equations 
Rw — 39 R= —KT yy 
with 
Lun = L (p Yury — Pry Vepr y Vr Pru Pru Vv Y). (37) 


1.5 The action with the y and Y fields not coupled 


In this section we shall use the action (cf. (31)) 
[{Tr (8° Op) — 2% (Gy* Pa Guy p)} dee (38) 


which contains derivatives up to the first order only of the quantities Y, P, vy" 
and Dy These will be varied independently, “ according to (30) and I’, under 
the restriction (34), i.e.2 : 


: The first part of this action was proposed by Klein in [1]. 
One could fear to lose, by this restriction, some equations obtained from variations of Tu 
corresponding to S transformations and needed for the antisymmetric part of the ‘“€,-equa- 
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4) ‘fot Py Tot think ‘edt: we write 
T= Pit My, tig ag (41) 


re My is a tensor-operator from which we require ones it fulfills 


Lage te clr a=), 1 (42) 


ge 


: pete ie Eicane 
[s, M]= — ty", pp + oP YY (43) 


oe M, can be solved from equations (42) and (43) by expanding it in the 7 
- matrices of Appendix 1 and determining the coefficients. The result is 


v 


whe es 


oN Te RR 


Be ie ie I a 


aN ON 


PRR eR 


Es 


( a =-1 = a 
My = {BP Sa VEPs tI PY 58a Ps YS PY Yul + Cus, (44) 


where the fact was used, that a three-index quantity symmetric in one pair and 
antisymmetric in another pair of indices must vanish. The c-number, weight —1, 


vector density, Cu, represents the only arbitrariness of the solution. It is convenient 


to give expression (44) without the last term a special symbol 
Myu=M ap Crys: (45) 


It is seen by inspection of (44) that 
[y", my] =9. (46) 


tions” (see eq. (30)) to vanish automatically. However, the only “§ variation” (see (A2.10)) 
equation lost by the restriction (39) is easily seen to be 


(PY* P),a = 


and this equation is obtained also from the variation of p and @ (see (57)). 
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pt : 


yiver ne +n gitcibaoes 2 


1 1 sy tikes ay 
hats ; , Pu 7up=0 aie porerteetr rey : 
and iW Puy +P y'u= 0. i, wil (Gb) evibe ot datw # 


Papedicine the new I’ field 


a Tt =~ v= Cuys 


we see, with the help of (47), that (49) and (50) ene 
 Pru=O 
and Pre y*=0, ° x : 


in which equations the arbitrary field c, is canceled. 
Equations (52) and (53) can also be written, using the ry field, as 


Yu My =0 
and Pru y “Ee y M= 0, 
where May hy", Mus =31y", my} = y" mp = My". 


Equations (45) and (46) show the validity of the relations appearing in (56). 
Equations (54) and (55) show how the new Dirac equations differ from the 
ordinary ones, in which the last terms containing the matrix M are not found. 
In Appendix 4 a relatively simple expression for MV is given. Equations (54) and 


(55) imply 
(Fy"g)u=0 (57) 
that is, the charge continuity equation. 


ce. Variation of y" 


Varying y“ and 8” in (38) according to (A2.1) and (A2.2) gives 


™ 


di (Diy s i 50uy ®,, 8°”) ae (p Va Yin Pu Va wy) oF (ot, m1) =0. 
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_ If we wish to compare these equations with Einstein’s equations it is conven- 
lent to change the e.c.d.-s with respect to I, to e.c.d.-s with respect to I’. 
This can be done without changing anything else in the equations. Furthermore 
_ we can use (48) in the first two terms noticing that the second term on the right- 
' hand side is antisymmetric. The result is 


: Ruy > $9Iuv k= 
cad 7 <r ame, —_ 
= GP Ye Pru Pru YeY — Pry Yup +P yu prrr) — 


te ies 
SPAY, Mod Yt 7 Iu P(e, hy. 


In this expression the second term on the right-hand side cancels exactly the 
contribution of m, to the first term (notice that the second term is symmetric 
according to (A3.3)), so that the e.c.d.-s with respect to I’ appearing there are 
changed to e.c.d.-s with respect to I’. We obtain finally 


R= 2G R= SRL 39m P My), (58) 


where 7',, and M are the quantities defined in (37) and (56) respectively. 
Equation (58) shows how the “new Einstein’s equations” differ from the ordi- 
nary ones, in which the term containing M does not appear. The new equations, 
(58), can be given a remarkably simple form if one notices that according to 
(37) and (54)—(55) 
T*,=ypMy. (59) 


With the-help of (59) the equations (58) obviously reduce to 
Big HH L yy (60) 


Hence, when x=0, the action principle employed here gives exactly Einstein’s 
equations in vacuum. 


II. Conservation theorems 


In [5], Moller has shown that his well-known conservation theorems for energy 
and momentum in the theory of general relativity (see Moller [9]), in the form 
where —x7',, has been replaced by the left-hand side of Einstein’s equations, 
are the identities derived from local varitions of the g,, in the curvature 


scalar. 
In this section we shall treat the gravitational part of the Lagrangian in (31), i.e. 


Tr (3° Op), (61) 


in a corresponding manner, by varying y” and I", locally (cf. Laurent [6)). The 
local variations of a few quantities are collected in Appendix 2. There and hence- 
forth the local variation is denoted by the symbol 0. 
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In this section no other I’ field than the I” field will be used. For simplicity : 
we therefore omit the prime connected with Ty, everywhere, i.c. 


are a0 (62) 


II.1. Conservation theorems connected with coordinate transformations. 


We shall here investigate how (61) changes under a local variation connected 
with the infinitesimal coordinate transformation 


alt = al + EM (x) (63) 


i.e. a variation of type (a) in Appendix 2. 
Using (10) and (62) followed by (A 2.6) and a few partial integrations we obtain 


6 | Tr 8" ®,,)d*x= 
= =f Tr {(O 8%) 9+ (Dug 3"2),p — (Ppp 3”) ,6 + Duy 3” o} Ed* 2e. 


This variation of (61) vanishes, however, due to the invariance of (61) under 
coordinate transformations and we obtain the identity 


Tr {(O.,3°)o—s Dien? 5 =, (64) 
This identity is closely related to the identity 
(Rx — j g% N).5 ==() 


which is obtained in a corresponding manner from the action i Rd*2=0. Fur- 
thermore 


6 Tr (®,, 3") = — Tr (By, 3” é2),, (65) 
due to the fact that (61) is a scalar density. On the other hand is 
6 Tr (®,, 3”) = 
= OTe {Dy e+ Dg £2) hy. + Te {Dy (8 Eg + 8" Eg — 8" Eg —8""n£°)}, (66) 


where (10), (A2.7) and (A2.6) have been used. 


_ Identifying the coefficients for £ in the right-hand sides of (65) and (66) 
gives the following relation 


ov = he B= 0, (67) 


where we have made use of (64). 
Identifying the coefficients for &, gives 


Tr(O5o" =e te ins (68) 


while the identification of the coefficients of £¢,, gives nothing. 
With help of the antisymmetric pases 
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APP TT BOM) i - (69) 
we can thus according to (68) and (67) define a divergence-free quantity 
TP= AP, (70) 
ag p= 0. (71) 


In the following section we shall discuss the relation of this quantity to the 
Moller energy-momentum complex. 


II. 2. Discussion of the new superpotential 


_The superpotential A,*” in (69) may be written in a different form 
A,P? = —4VgTr (ty", Pel’) = 
= WgTr(27")= 
= eG) F209 


The superpotential A,*” is, however, antisymmetric in fb and y. It may thus 
be written 


i tt alles 9 2S Valeur adle g (72) 
where Hal = — Moa GY + Mn 9 (78) 


is Moller’s [9] superpotential. 

Our derivation of &,° is purely formal in nature and we have, of course, no 
right to call it the energy-momentum complex. One circumstance that casts doubt 
on such an interpretation is that neither ZT,’ nor A,” is invariant under S 
transformations. 

It is tempting to guess that the last term in (72) can be completely trans- 
formed away by an S transformation, in which case A,” would be ‘‘almost”’ 
identical with Moller’s superpotential. To investigate if this is the case seems, 
however, to be a rather complicated undertaking. Let us instead briefly discuss 
the properties of A,’” which can be seen immediately. 

From (69) it is clear that A,’’, in distinction to 7,°’, is a tensor density and 
from (72) it follows that Ao®” is invariant under time-independent S transfor- 
mations. It is also clear from the same formula that Ao” goes over into Moller’s 
Yo” in the case of a static y field. This in turn has as a consequence that 
Zp° integrated over the space of a closed system gives the correct result. 

It is also interesting to note that the essential property of the Moller 2o° 
that it transforms as a vector density under purely spatial transformations, is 
a property of Ao’, also. Equation (72) shows this immediately. 

It is perhaps not necessary to point out that A,’”,, transforms as a tensor 
density under linear transformations. 


ry a 
20 > 
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Omer ninse brermafr wei rite 


O Be - Pail j? 
“s gives, with, tho mee of oy: oa (A29), cm ipl 4 ee Pia : an © igs: iL ® 
ww sod qaiwo ies oft arr 
_ [®,», 3”]=0. nye adie rye VyEure wie? 
—-‘Tt further follows from _ = — 
6 Tr (O,, 3”) =0 
ais 
: that, following the same pattern of derivation as in section II.1, but using for- 
mulae (10), (A2.9), (A2.10) and (74), we obtain gt 
(3, T.]|4=0 (75) 
; pe 8 se are (76) 
: : “{ : 
; Conclusion 


The purpose of the first section of this paper was to investigate the changes. 
in the Dirac and Einstein equations when one breaks the relation (12) between 
the y and I’ matrices in the combined action (31) for the gravitation and Dirac 
fields, and varies the two sets of matrices independently. The result are contained 
in formulae (54), (55) and (58). It should be pointed out that these equations 
differ from the Dirac and Einstein equations only by small terms of the order 
x* and any experimental verification of one set or the other is at present far 
beyond reach. Now when it is clear that the variational principle leads to a 
reasonable result, the main interest would probably lie not so much in the second 
order equations as in the action itself, which may be used subsequently as the 
starting point for a quantization procedure. 

In the second section we showed the simplicity of the procedure for obtaining 
conservation identities when one starts with the gravitational part of the action 
(31). In particular the superpotentials obtained are automatically antisymmetric. 
We also saw there that the new superpotentials are closely related to those 
introduced by Moller. 
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Tr (y"y")= 49", 
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bey = 29. 
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(A1.7) 
(A1.8) 


: ; ey: i . 
Pantha cgay aL Se, enstaye tan 
.. + ee ra ~— ae ae 
tidy! @ygt ce mlesmt % weit 
= gives, for the variation $=? 2”,,.6 = 1" eP 
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O38" =\C', ae nae - Sy Ca 


% 
The corresponding variation of the affine connexion is, as is well known, 


) Is, = 39 [(6 Gop): + (OGoy).8 — (6 Jpy):01 ; ae | : aid 


or, introducing the €",, q ~ 


cna 
-, 


OD, = — 9° [Ecp:y + Copia — Esy:o]- : (A 2.3) 
Using this formula we obtain 
OT, =4 Go. 8°? + 6a,. (A 2.4) 


Equation (A 2.4) follows, e.g., by comparing the expressions for 6 (y",,) and (dy), 
the first of which equals zero. The quantity a, is the vector field appearing in (13). 
Below we have collected local variations corresponding in (a) to the infinitesimal 


coordinate transformation w“=a"-+ &"(a), and in (6) to the infinitesimal S trans- 
formation with S=1+o. 


(a) Cpe Oe rae (A2.5) 
53M Sg EF a gre ce gue a gee (A 2.6) 
Ole + yay ere (A 2.7) 
(b) éy"=[y",.0] (A2.8) 
53” = [8", o] (A 2.9) 
éT,=(C ole (A 2.10) 
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_Tr(e"[M,, 2 NY oat 
= —Tr ([s°’, M,] M,) 
=xTr ({*, pp} ML) = 
- =x ptM,, y*} y. 


~ ~ Rn 2 gue 


pe #4 
Hence the complete expression reads 
hoe . wn 


Tr (98° Dep) = 2B + (PL8% 3} at 2B ltn 7h Ws (A3.2) 


where the comma in the middle term is changed to a semicolon so that (A3.2) 
becomes covariant. Further M, has been changed to m, in the last term. We are 
- allowed to make this change, according to (45) and (A 1.6). 
: We notice that the last term in (A3.2) is symmetric in y and @g, so that 
J 
e 


FE Pim, Yu p= Plime Yo} Y- (A3.3) 
This follows immediately from (44). 


APPENDIX 4 


In this appendix we shall simplify the explicit expression for M, defined in 
(56). Using (44) in (56) gives 


M== bP {5x0, vy y (8, v3 — 29 Ps Sap vss —$ Py (A4.1) 


i 
8 
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